Monte Carlo simulations of the three-dimensional XY spin glass focusing on the chiral 

and the spin order 



O 

(N 

o 

Q 



C/3 



o 
o 



o 

oo 
m 



Tomoyuki Obuchi and Hikaru Kawamura 

Department of Earth and Space Science, Graduate School of Sctence, 
Osaka University, Toyonaka, Osaka, 560-0043, Japan, 
(Dated: December 18, 2012) 

The ordering of the three-dimensional isotropic XY spin glass with the nearest-neighbor random 
Gaussian coupling is studied by extensive Monte Carlo simulations. To investigate the ordering of 
the spin and the chirality, we compute several independent physical quantities including the glass 
order parameter, the Binder parameter, the correlation-length ratio, the overlap distribution and the 
non-self-averageness parameter, etc, for both the spin-glass (SG) and the chiral-glass (CG) degrees 
of freedom. Evidence of the spin-chirality decoupling, i.e., the CG and the SG order occurring at 
two separated temperatures, Tea > Tsg > 0, is obtained from the glass order parameter, which is 
fully corroborated by the Binder parameter. By contrast, the CG correlation- length ratio yields a 
rather pathological and inconsistent result in the range of sizes we studied, which may originate from 
the finite-size effect associated with a significant deviation of the spatial CG correlations from the 
standard Ornstein-Zernike form. Finite-size-scaling analysis yields the CG exponents vcG = 1.4±0.1 
and ricG = 0.29±0.12, and the SG exponents vsa = 1.23to ll and rjsG = -0.42j;E;;^^. The obtained 
exponents are close to those of the Heisenberg SG, but are largely different from those of the Ising SG. 
The chiral overlap distribution and the chiral Binder parameter exhibit the feature of a continuous 
one-step replica-symmetry breaking (IRSB), consistently with the previous reports. Such a IRSB 
feature is again in common with that of the Heisenberg SG, but is different from the Ising one, which 
may be the cause of the difference in the CG critical properties from the Ising SG ones despite of a 
common Z2 symmetry. 



I. INTRODUCTION 

In spite of long history of research, spin glass (SG) is 
still a hot topic in statistical physicsi. It is a typical 
system possessing both strong frustration and random- 
ness, leading to several extraordinary behaviors such as 
the slow dynamics and the rejuvenation-memory effect. 
In a theoretical treatment of SG, Edwards and Ander- 
son (EA) proposed as early as in 1972 a simple modeP, 
the so-called EA model, in which the spins are put on 
each site of a regular lattice and interact via the random 
coupling taking both positive and negative signs. The 
infinite-range or the mean-field version of the EA model, 
first presented by Sherrington and Kirkpatrick (SK)2,, 
was solved by Parisi revealing an intriguing concept of 
the replica symmetry breaking (RSB)'*. For both cases 
of the Ising and the Heisenberg SK models, the relevant 
RSB turned out to be of hierarchical nature. In spite of 
such success of the mean-field theory, understanding the 
nature of the ordering of the finite-range EA model in 
three dimensions (3D) still remains incomplete. Numeri- 
cal simulations have been the main tool in attacking the 
issue. Although the existence of a finite-temperature SG 
transition was established in the 3D Ising EA model^"— , 
earlier numerical simulations on the 3D XY and the 
Heisenberg EA models suggested the absence of a finite- 
temperature transitiorii^^— , in apparent contrast to ex- 
periment^i^Ti^. 

Some time ago, in discussing the ordering of frustrated 
vector spin systems, Villain proposed a possible signifi- 
cance of the "chirality" degree of freedom, an Ising-like 
scalar quantity which represents the handedness of the 



noncoUinear spin structure. Villain made a conjecture 
that the 3D XY SG might exhibit a finite-temperature 
SG ordering by noting the Ising nature of the chirality 
and by invoking the occurrence of a finite-temperature 
SG transition in the 3D Ising SG^^. 

Numerical simulations on the vector SG, i.e., the two- 
component XY SG or the three-component Heisenberg 
SG, investigating both the spin and the chirality degrees 
of freedom, have been performed ever since 198522.. A 
crucially important concept which emerged from these 
studies is the possible "spin-chirality decoupling" phe- 
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. In 3D, it means that the chirality orders 
at a temperature higher than the spin, with an interme- 
diate "chiral-glass" (CG) phase where only the chiral- 
ity exhibits a glassy long-range order while the standard 
SG order still remains short-ranged. In dimensions lower 
than three, where both the spin and the chirality order 
only at r = 0, the spin-chirality decoupling means that 
the spin and the chiral correlation-length exponents are 
mutually different, i.e., the existence of two different di- 
verging length scales at the T = transition'^-"—. The 
concept of the "spin-chiraity decoupling" in SG was orig- 
inally proposed by one of the present authors (H.K.), and 
subsequently explored further by the corroborators for a 
wide range of systems. In terms of a symmetry, in the CG 
phase the Z2 spin-refiection symmetry is spontaneously 
broken with keeping the SO{3) or SO{2) spin-rotation 
symmetry unbroken. 

In case of the 3D XY SG, which is a target of the 
present study, the spin-chirality decoupling was first ex- 
amined by the numerical domain-wall renormalization- 
group calculation and also by a Monte Carlo (MC) sim- 
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ulatio n'^^' 1 . Further interesting features revealed by 
these numerical analyses might be the possible one-step 
RSB (IRSB)-like feature of the CG ordered state, and 
the non-Ising character of the CG criticality^. 

Similar spin-chirality decoupling phenomena, the 
IRSB-like nature of the CG ordered state and the non- 
Ising character of the CG criticality were also observed 
in the 3D Heisenberg SG, a reference model of many of 
realistic SG materials including canonical SG, in spite 
of the difference in the nature of the chiralities relevant 
to the XY and the Heisenberg spins. For the former, 
the chirality is quadratic in spins and time-reversal even, 
while, in the latter, it is cubic in spins and time-reversal 
odd. Indeed, on the basis of the spin-chirality decoupling 
picture of the 3D Heisenberg SG, the chirality scenario 
of the experimental SG order was advanced22r^. Recent 
large-scale simulations have revealed that the SG order 
actually takes place at a nonzero temperature^^—, in 
contrast to earlier beliefs. Besides, most of recent simula- 
tions point to the occurrence of the spin-chirality decou- 
pling in the system^"—, i.e., < Tsg < Tea- Several 
experimental facts were also successfully explained by the 
chirality scenario^SiSi, and it becomes increasingly clear 
now that the spin-chirality decoupling is an indispensable 
concept in understanding the realistic SG systems. 

On the other hand, the situation in the 3D XY SG 
seems less clear. Extensive calculations comparable in 
their scale to those of the 3D Heisenberg SG have been 
scarce, and the occurrence of the spin-chirality decou- 
pling still remains controversial. Maucourt and Grem- 
pel suggested on the basis of their T = domain- 
wall renormalization-group calculation for lattices L < 8 
the occurrence of a nonzero Tsg located below Tcg^- 
Mentioning some of the recent MC simulations on the 
model, Kawamura and Li simulated the ±J EA model 
by an equilibrium MC simulation up to the linear size 
L = 16, suggesting the occurrence of the spin-chirality 
decouplings^, whereas Granato performed the dynam- 
ical Langevin simulations of the model for lattices of 
L < 12, to conclude the occurrence of a single transition 
Tsg — Tgg^'^- Nakamura and collaborators performed 
the nonequilibrium relaxation analysis for lattices up to 
L — 55 (this method enables one to treat relatively larger 
sizes but some drawbacks appear in its short-time obser- 
vations), and suggested that Tsg and Tcg were identical 
or close even if they were to be different^i^. Young and 
corroborators investigated the Gaussian EA model by 
equilibrium MC simulations for lattices up to L < 24, re- 
porting no evidence of the spin-chirality decouplingS2i^. 

This confusing situation motivates us to re-examine 
the ordering of the 3D XY EA model with the ran- 
dom Gaussian coupling by large-scale MC simulations 
by treating large sizes up to L = 40, considerably larger 
than the sizes studied before by equilibrium simulations. 
A large number of samples of order Ng ^ 0(10"^) are 
simulated to obtain reasonable statistics. Furthermore, 
we compute various independent quantities including the 
glass order parameter, the Binder parameter, the corre- 



lation length ratio, the overlap distribution function and 
the non-self-averageness parameters both for the SG and 
the CG, in order to check consistency among various in- 
dependent quantities. 

We note that the 3D XY EA model is a reference 
model for SG magnets with an easy-plane-type uniaxial 
magnetic anisotropj*^"— . Readers are referred to Ref. 
74 for detailed discussion. The ordering properties of the 
model would also be helpful in understanding the pecu- 
liar ordering behaviors experimentally observed in these 
granular cuprate superconductors'^^—. 

Overall, the results of our large-scale simulations speak 
for the occurrence of the spin-chirality decoupling in the 
3D XY SG. The estimated SG and CG transition tem- 
peratures are T, = 0.2741° °^^ ^^^^ ^ o.308 ± 0.05, 
Tcg being higher than Tsg by about 10%. In estimating 
the transition temperatures, we have found a reasonably 
good consistency among various independent quantities, 
with one exception of the CG correlation-length ratio, 
which behaves rather badly leading to a pathological es- 
timate of Tcg- Thus, in deriving the above estimate of 
TcGi we have not used the CG-correlation length data, in 
contrast to Ref. 67. To clarify the origin of the observed 
pathological behavior of the CG correlation length, we di- 
rectly compute the spatial chiral correlation function to 
find that the standard Ornstein-Zernike (OZ) form may 
be inappropriate to describe the CG correlation length in 
the range of small sizes in which we make our simulations. 

The critical properties associated with the SG and the 
CG orderings are also examined. We obtain the CG ex- 
ponents i/cG = 1.4±0.1 and rjcG = 0.29±0.12, while the 
SG ones vsg = 1 -2310:^6 and rjsG = -0.421°;^?, where 
V and r] are the correlation-length and the critical-point- 
decay exponents, respectively. These exponents turn out 
to be close to the corresponding Heisenberg SG expo- 
nents, but are different from the Ising SG ones. We also 
confirm the IRSB nature of the ordered state, consis- 
tently with the previous reports. 

This paper is organized as follows. In Sec. [Hi we in- 
troduce the model and explain some of the details of 
our simulation. In Sec. Illli we define several physical 
quantities which we compute to examine the SG and 
the CG orderings. The results of our MC simulations 
are presented in Sec. IIVI We show the data of the glass 
order parameter, the Binder parameter, the correlation 
length ratio, the overlap distribution function, the non- 
self-averageness parameters and the spatial correlation 
function for both the SG and the CG degrees of free- 
dom. The CG and the SG transition temperatures are 
estimated via an infinite-size extrapolation of appropri- 
ate finite-size data. Possible RSB character of the or- 
dered state is examined in this section. In Sec. |Vl the 
critical properties of the SG and the CG transitions are 
analyzed, and the CG and the SG critical exponents are 
determined. Comparison is made with the corresponding 
exponents of the 3D Heisenberg and of the 3D Ising SGs. 
The last section is devoted to summary and discussion. 
In appendix, the behavior of the SG Binder parameter 
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in the thermodynamic limit across the CG and the SG 
transition points is analyzed. 

II. MODEL AND SIMULATIONS 

The model we study is the isotropic XYFiA model on 
a 3D simple-cubic lattice. The sites are labeled by the 
index i [i — 1, 2, • • • , N), the corresponding coordinate 
being denoted as = (xi,yi, Zi). The total number of 
spins N is related to the linear system size L a.s N = L^. 
The XY spin on the ith site, Si, has two components, 
St = {Six,Siy) = (cos 6*1, sin 6*^) where < 9i < 2tt. The 
Hamiltonian is given by 

H = — JijSi ■ Sj, (1) 

where the summation (i, j) is taken over the all nearest- 
neighbor pairs. The interaction Jij is a random Gaussian 
variable whose mean and variance are taken to be zero 
and unity, respectively. The partition function is given 
by 

where /? is the inverse temperature 1/T normalized by 
the Boltzmann constant fee ■ The thermal average will be 
denoted by the angular brackets (•••). 

We perform MC simulations based on the single- 
spin-flip Metropolis method combined with the over- 
relaxation method and the temperature-exchange tech- 
nique. This algorithm is known to effectively reduce 
the long correlation time involved in simulation of hard- 
relaxing systems such as SG. 

In a unit process of the over-relaxation, we compute 
first the local field felt by the spin at site i, hi — 
X^jeA '^ijSj where represents the neighbors of the 
site i, and then refiect the spin Si with respect to the 
local field hi as 

S.,^S[ = -S, + 2^^h,. (3) 

The simple cubic lattice consists of two interpenetrat- 
ing sublattices, and we perform the Metropolis update 
sequentially through the sites on one sublattice after an- 
other, which is followed by the M-times over-relaxation 
sweeps, also performed sequentially through the sites on 
each sublattice. This procedure consists our unit MC 
step. In our simulations, we take M equal to the linear 
system size L. 

In the temperature-exchange process, we prepare Nt 
spin configurations at a set of temperatures distributed 
between Tmin and Tmax- The maximum temperature 
Tmax is chosen to be high enough where the autocorre- 
lation times of the spin and the chirality are sufficiently 
short even in the single-spin-flip dynamics, typically 40 



MC steps per spin (MCS) . Each trial of the temperature 
exchange is performed for a pair of neighboring temper- 
atures. A temperature-exchange trial is done after every 
MCS. 

In Table HI we summarize the simulation parameters 
employed in our MC simulations, which include the lin- 
ear system size L, the total number of samples aver- 
aged (independent bond realizations) iV^, the number of 
MCS discarded for equilibration Nmci, the number of 
MCS employed for measuring physical quantities Nmc2, 
the maximum and the minimum temperatures in the 
temperature-exchange process Tmax and Tmin, and the 
number of temperature points Nt- In simulating SG 
systems, special care has to be taken for equilibration. 
In this paper, we follow the criteria of Ref. 55. Error 
bars are estimated by using the bootstrap method from 
sample to sample fluctuations. 



L 


Ns 


Nmci 


NmC2 


7m ax 


Tin in 


Nt 


4 


5000 


1 X 10" 


1 X 10" 


0.86 


0.24 


32 


6 


5000 


3 X 10* 


1 X 10" 


0.86 


0.24 


32 


8 


2000 


5 X 10* 


1 X 10" 


0.80 


0.24 


32 


12 


2000 


1 X 10^ 


1 X 10" 


0.60 


0.24 


32 


16 


2048 


4 X 10^ 


4 X 10" 


0.52 


0.26 


32 


20 


1024 


5 X 10^ 


5 X 10" 


0.50 


0.266 


40 


24 


1024 


7.5 X lO"^ 


7.5 X 10" 


0.49 


0.271 


40 


32 


1024 


1.5 X lO** 


1.5 X 10*^ 


0.48 


0.2736 


56 


40 


384 


2-2.8 X lO** 


2-2.8 X 10*^ 


0.46 


0.2891 


64 




128 


2-3.4 X 10® 


2-3.4 X 10*^ 


0.442 


0.2792 


64 



TABLE I. Parameters of our Monte Carlo simulations. L 
is the linear system size, Na is the total number of sam- 
ples, Nmci is the Monte Carlo steps per spin discarded for 
equilibration, Nmc2 is the Monte Carlo steps per spin subse- 
quently used in measuring physical quantities, Tmax and Tmin 
are the highest and the lowest temperatures employed in the 
temperature-exchange process, and Nt is the total number of 
temperature points. For L = 40, we use two different temper- 
ature sets, where Nmci and Nmc2 are adaptively chosen for 
each set to satisfy the equilibration condition. 



III. PHYSICAL QUANTITIES 

In SGs, the conventional order parameter is an overlap 
between two independent systems with a common Hamil- 
tonian. In the case of the XY model, each spin has two 
components and the spin overlap becomes a tensor with 
indices a and /3 (a,/? = x,y). We deflne the wavevector 
fc-dependent spin overlap as 

1 ^ 

9"Mfc) = ]^E^-^l?e*-, (4) 

i=l 

where the superscripts (1) and (2) denote two indepen- 
dent systems with the same Hamiltonian. For simplicity 
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of notation, we write 



Qs{k) = ka/3(fc)r 



(5) 



Similarly, we define the chirality and introduce the as- 
sociated overlap. The chirality at a plaquette p which is 
perpendicular to the /i(= x,y,z) axis is defined by 



— y 



sgn(Jy )sin(6'i - 6*^), 



(6) 



where the directed sum V!/- is taken over four 
bonds surrounding the plaquette p in a clockwise direc- 



tion- 
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I. The chiral overlap is then given by 

^« ~ iv ^ p-L^ 



(7) 



The SG order parameter qgQ and the SG susceptibility 
XSG are then defined by 



(8) 



where the square brackets [• • •] denote the configurational 
average, i.e., the average over the bond disorder. 

The SG Binder parameter gsc and the SG correlation 
length S^sG are defined by 



gsG 



isG = 



1 [{qsW)] 

2sin(fc„,i„/2)V [((Z.(fc,ni„)2)] 



-1, 



(9) 
(10) 



where fc,nin = {2tt/L, 0, 0). Note that the OZ form of the 
correlation is implicit in this definition of the correlation 
length. 

On the other hand, the CG order parameter and the 
CG susceptibility are given by 



<i-ci = Mm] , xcG 



CG- 



(11) 



The direction /i-dependence of the right-hand side should 
vanish after the sample average [•••], and we take the 
average over ^ = x,y, z in the actual calculation. The 
CG Binder parameter and the CG correlation length are 
given by 



1 



gcG 



■=GG 



_^ Mm] 

2sin(A:„i„/2)V [(|'?^^(femi„)P)] 



- 1. 



(12) 



(13) 



Although the /i-dependence again vanishes for gcd it 
remains for due to the nontrivial wavevector depen- 
dence of 9j^ (fcmin)i j-e., the dependence on the direction 
fi with respect to fc(|| x) taken here parallel with x. We 



denote as S^fjQ and as £,qq, and will show both 
the data below. 

We also define a parameter quantifying the non-self- 
averageness of the order parameter, the A parameter—. 
It is defined either for the spin or for the chirality by 



A 



SG 



AcG = 



mow 



(14) 
(15) 



The /Lt-dependence vanishes for Acg- The A parameter 
becomes nonzero if the SG or the CG susceptibility is 
non-self-averaging. Note that, even when qsG vanishes 
in the thermodynamic limit, AgQ can become finite if 
XSG is non-self-averaging. In the current problem, such 
a situation can emerge in the temperature range Tsg < 
T < TcG in the possible occurrence of the spin-chirality 
decoupling. 

We also introduce the so-called Guerra parameter or 
the G parameter defined by 



GsG — 



GcG — 



{qsiorY -[{qsmi' 



[(g.(o)4)] - [{q4or)Y 



mom - mor)r 



(16) 
(17) 



The G parameter looks like the A parameter, but there 
is a difference in that the G parameter can be finite 
even when the ordered state does not accompany the 

RSB84-86. 

The distributions of the spin and the chiral overlaps 
might provide a signal of the RSB. In this paper, we 
examine the following two overlap distributions 



Ps{q)= <5(^g-^9c.o(0) 

p.{q) = [Hq-qm)]- 



(18) 
(19) 



For the chiral overlap distribution P^iq), the RSB ef- 
fect is simple: if there is no RSB in the ordered state, 
the distribution has only two (5-peaks in the thermody- 
namic limit which are related each other by the Z2 spin- 
reflection symmetry of the whole spins. On the other 
hand, the spin overlap distribution Psiq) takes a non- 
trivial form even in the SG ordered state without the 
RSB, a superposition of two (5-peaks located at g = iqEA 
and a broad distribution spanning between these diverg- 
ing peaks, due to the projection of the tensor onto 
the diagonal component. This makes obtaining a clear 
indication of the RSB from the Pg data rather difficult. 
For further details, see Ref. 47. 
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IV. MONTE CARLO RESULTS 

In this section, we present the resuh of our MC simula- 
tions. We show first the temperature dependence of the 
specific heat in Fig. [TJ No appreciable anomaly is seen in 
the specific heat, though the SG and the CG transition 
points actually exist in this temperature range, as will be 
shown below. The CG and the SG critical temperatures 
TcG and TsG are denoted by arrows in the figure. 




I ' ' 1 ' ' ■ ' ' ■ — 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

T 



FIG. 1. (Color online) The temperature and size dependence 
of the specific heat per spin. Magnified view is given in the 
inset. The arrows at higher and lower temperatures indicate 
the CG and the SG transition points, respectively. 

To investigate the SG and the CG orderings, we show 
the temperature dependence of the SG and the CG sus- 
ceptibilities, XSG and XCG, in Fig. [5] In contrast to the 
SG susceptibility xsG, which tends to increase as the sys- 
tem size L is increased in the whole temperature region, 
the CG susceptibility XCG shows such a behavior only in 
the temperature region T ^ 0.4, whereas, in the region 
T > 0.4, it exhibits an opposite size dependence. This 
implies that the CG critical region is relatively narrow, 
which is common with the observation of the 3D Heisen- 
berg SG^. 

In Fig. [21 we plot the size dependence of the SG and 
the CG order parameters q'gQ and q^^]-, for several tem- 
peratures on a double-logarithmic plot. The data of 
the CG order parameter q)jQ exhibit a straight-line be- 
havior around a temperature T ~ 0.306. It exhibits 
a clear upward trend at lower temperatures implying 
the appearance of the CG long-range order, whereas at 
higher temperatures it exhibits a downward trend, even- 
tually approaching another straight line with a slope 
—d = —3 generally expected in the disordered phase for 
large enough systems. Then, we get our first estimate 
of the CG transition temperature, Tcg ^ 0.31. By con- 

(2) 

trast, q^Q exhibits such an upward trend only at the low- 
est temperature studied, T = 0.266, with a straight-line 
behavior observed around T ^ 0.276. Then, we get an 
estimate of the SG transition temperature, Tsg ~ 0.28. 



Hence, our data of the size dependence of the glass order 
parameters g*-^-* suggest that the spin-chirality decoupling 
really occurs in the present model. 

In Fig. m we show the SG and the CG Binder pa- 
rameters. Consistently with the earlier reports**^'^^, the 
CG Binder parameter exhibits a non-divergent dip and 
a crossing among different sizes on the negative side of 
gcG- Such a behavior is expected in a system exhibit- 
ing a continuous one-step RSB (IRSB). The crossing and 
the dip temperatures are expected to converge to Tgg in 
the thermodynamic limit, which might provide a way to 
precisely estimate Tgg- By contrast, the SG Binder pa- 
rameter gsG shows no crossing nor dip, decreasing mono- 
tonically as the system size is increased. So far, informa- 
tion concerning the transition points has been hard to 
obtain from gsG- We re-examine the relation of gsG to 
the CG and the SG transition points below in appendix, 
and point out a possibility to extract information about 
Tsg and Tcg from the gsG data. 

In Fig. O we show the temperature dependence of the 
SG correlation- length ratio (,sg/L. Those of the parallel 
CG correlation-length ratio ^qq/L and of the perpendic- 
ular CG correlation-length ratio £,c'g/^ ^^'^ given in Fig. 
[6l Both the SG and the CG correlation-length ratios 
show crossing among different sizes, and the crossing tem- 
peratures are expected to converge to the corresponding 
critical temperatures in the L — >■ oo limit. 

To estimate Tsg, we plot in the left panel of Fig. [7] the 
crossing temperatures of the SG correlation-length ratio 
^sg/L for pairs of sizes L and sL with s = 2, 3/2, 4/3, 5/3 
and 5/4 versus 1/Lave where Lave = (1 + s)L/2. Note 
that the crossing temperature for the pair (L, sL) = 
(32,40) with s = 5/4 is estimated by extrapolating the 
data to lower temperatures, since the raw data of S,sg/ L 
do not show a crossing in the investigated temperature 
range. 

We then try an infinite-size extrapolation of the cross- 
ing temperature Tcross{L; s) based on the scaling from, 

T„ossiL;s) = T, + CsL-'', d = ^ + \- (20) 

In the fit of Tcross{L; s) of the SG correlation-length ra- 
tios S^sg/L, we perform a combined fit of different s- 
sequences with a common Tc = Tsg and a common 
e = OsG- We use in the fit the data for s = 2, 3/2, 4/3 and 
5/3 only, not the one for s = 5/4 because of a possible in- 
accuracy due to the extrapolation employed for the data 
of (L, sL) — (32, 40) mentioned above. The resulting fit- 
ting curves are shown in the left panel of Fig. [T] The 
optimal fit is obtained for Tsg — 0.274 and Osg — l-l- 
To estimate the associated error bar, we plot in the in- 
set of Fig. [8] the total x-square value of the fit versus the 
assumed Tsg value. The horizontal line in the figure rep- 
resents the total x-square value greater than the optimal 
value observed at Tsg = 0.274 by unity, which gives our 
error criterion. The asymmetry of the curve results in the 
different values for the upper and the lower error values. 
We then quote Tsg = 0.2741°;°^^^ and Osg = l-lOloJs- 
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FIG. 2. (Color online) The temperature and size dependence of the SG susceptibility XSG (left), and of the CG susceptibility 
XCG (right). The insets are magnified views. As can be seen from the insets, the magnitude of xsG increases as the system 
size is increased in the whole temperature region, while that of xcG exhibits an opposite size dependence in the temperature 
region T > 0.4. 




58 



10 12 16 20 24 

L 



40 48 



0.01 

0.004 
0.002 
0.001 

0.0004 
0.0002 
0.0001 




T=0.266 
T=0.276"--»' 
T=0.286---« 

T=0.296' B 

T=0.306 • — 
T=0.316— * 
T=0.326'- — 
T=0.336— 
T=0.346-- 



slope=-3 ■•. 



10 12 16 20 24 32 40 48 



FIG. 3. (Color online) The size dependence of the SG order parameter (left), and of the CG order parameter (right) on a 
log-log plot for several temperatures. Straight lines are drawn by fitting the three data points of smaller sizes of L = 4, 6, 8. 
The L = 24, 32 data for T = 0.266 as well as the L — 40 data for T — 0.266 and 0.276 are obtained by extrapolating the higher 
temperature data to lower temperatures. A line with a slope — d — —3, which is the expected large-L asymptotic behavior in 
the disordered phase, is also drawn. 



Next we turn to the estimate of Tcg based on the 
CG correlation-length ratio £,cg/L. Although the par- 
allel one ^qq/L and the perpendicular one ^qg/-^ S^^^ 
two different sequences of the crossing temperatures, they 
tend to accord for Lave ^ 8 a.s can be seen in the right 
panel of Fig. [7] (see also the inset). Hence, to extract 
Tcg, we use the data of for Lave > 8 only. Unfortu- 
nately, and somewhat unexpectedly, the result of the fit 
of £,cg/L turns out to be rather pathological. The esti- 
mated best TcG-value becomes extremely small or even 
negative. The resultant fitting curves are shown in the 
right panel of Fig. [7] (the lower curves) . We also exam- 
ine the possible change in the fit with varying the lowest 



size used in the fit, but the pathology cannot be cured. 
The cause of this pathology is not entirely clear, but may 
be due to the peculiar behavior of the chiral-glass corre- 
lations. For example, xcG shows a non-monotonic size 
dependence absent in the corresponding xsG- To further 
clarify this point, we discuss the properties of the spatial 
CG correlation function below in this section. Anyway, 
in the present analysis, we abandon the crossing temper- 
atures of ^cc/i in estimating Tcg- 

In estimating Tec, the crossing and the dip tempera- 
tures of gcG can also be used, and the data are plotted 
in the right panel of Fig. [71 The crossing temperature of 
gcG obeys the scaling form Eq. (1^ with Tc = Tcg and 
fi* = (^CG = ^CG + ^I^CG, whereas the dip temperature 
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FIG. 5. (Color online) The temperature and size dependence 
of the SG correlation-length ratio £,sg/L. The inset exhibits 
the wider temperature range. 



is expected to follow the scalmg form, 

Td^p{L) = TcG + ciL-^ + C2L-''^°. (21) 

Conventionally, the sub-leading correction-to-scaling 
term, C2£^^^^, is dropped because it gives a relatively 
smaller contribution than that of the leading correction- 
to-scaling term ciL~'^/'^'^'^ for larger L. In the range of 
system sizes of our simulations, however, we need to take 
into account this correction term for describing the size 
dependence of Tcross^ since the leading term ciL~^^'^'^'^ 
should describe the increase of the dip temperature with 
the system size, which is not observed in our simulation. 
(Such an increase of Tdip{L) was indeed observed in a 
recent simulation of the 4D Heisenberg SG*— .) The be- 
havior originates from the fact that both the dip and 



the crossing temperatures of gcG should converge to a 
common value, TqCi each with an exponent 1/v and 
9{> 1/v), while the crossing temperature always lies 
above the dip one. 

The combined fit of the crossing and the dip tem- 
peratures of gcG based on Eqs. (I20l2ip with a com- 
mon TcG and a common 9cg yields Tcg — 0.308 and 
(^CG — 0.88. In our estimate of error bars, we again 
calculate the total ^-square value against the assumed 
Tcg value, and the result is given in the main panel of 
Fig. [8l Our estimate of the CG transition temperature 
is then Tcg = 0.308 ± 0.005 and that of the exponent is 
OcG = 0.88 ± 0.03. Note that the leading term in Eq. 
(PT|) tends to be masked by the correction term, i.e., ci 
tends to be considerably smaller than C2, and the value 
of vcG is less precise. 

The SG and the CG transition temperatures Tgc = 
0.274t°:gJ2 and Tec = 0.308 ± 0.005 estimated from the 
Binder parameter are well consistent with the values ob- 
tained from the order parameter given in Fig. [31 Our 
error analysis indicates that Tsc and Tec are indeed 
separate. These observations certainly speak for the oc- 
currence of the spin-chirality decoupling in the 3D XY 
SG. The difference in the CG and the SG transition tem- 
peratures is about 10 percent, which is comparable with 
the corresponding value of the 3D Heisenberg SG*^^. 

In Fig. [21 we plot the ratio of the CG and the SG 
correlation lengths, £,cg/£,sg- For smaller sizes of 6 < 
L < 12, the ratio curves are almost size-independent^^ at 
lower temperatures, while for intermediate sizes of 16 < 
i < 24 the curves start to splay out but the tendency 
is still small^. For larger sizes of L = 32 and 40, the 
tendency becomes stronger and the ratio exceeds unity 
at low temperatures, which supports the spin-chirality 
decoupling ansatz. The intersection point of the ratio 
curves between different sizes lie around T ~ 0.31, which 
is consistent with our estimate of Tec = 0.308 ± 0.05. 
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FIG. 7. (Color online) The crossing temperatures Tcross of the SG correlation-length ratio ^sg/L of the two sizes L and sL 
are plotted versus the inverse mean lattice size 1/Lave where Lave = (1 + s)L/2 (left). The crossing temperatures of the CG 
perpendicular correlation-length ratio ^ca/^ ^^'^ of the CG Binder parameter gca as well as the dip temperature Tdip of gca 
are plotted versus 1/Lave (or 1/L) (right). The insets exhibit the wider size range. 



Now, to get further insight into the cause of the patho- 
logical behavior we encountered in ^cG: we discuss the 
CG and the SG spatial correlations, particularly how 
they obey to or deviate from the standard OZ form. The 
SG and the CG spatial correlation fvmctions, Cs{x) and 
Ci^ix), are defined by 



(22) 
(23) 



The computed Cg and for i = 16 lattices under pe- 
riodic BC are shown in Fig. [TU] for several temperatures. 
Note that the CG correlation function is normalized 
by its local amplitude so as to give unity at x = 0. The 



data are plotted on a scmi-logarithmic scale so that the 
standard OZ form gives a straight line. The leveling-off 
of the data observed at larger a; is a finite-size effect due 
to the imposed periodic BC. One can see from the fig- 
ure that Ck(x) drops rapidly from unity in the small- a; 
region of a few lattice spacings, strongly deviating from 
the OZ form, and becomes by an order of magnitude 
smaller than Cs{x) for larger x even below Tcg- Such a 
sharp drop of the spatial correlations at short length is 
much less pronounced in Cs{x). This feature of the CG 
correlations might cause some problems in defining the 
finite-size correlation length S^cG based on the OZ form, 
at least for small lattices treated in our present simula- 
tion. This is because the OZ-based definition of CcG, Eq. 
P^. contains the fc = part [((7k(0)^)], which is essen- 
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FIG. 8. (Color online) The total x-square values associated 
with the combined fit of the crossing and dip temperatures 
of the CG binder parameters gcG are plotted versus the CG 
transition temperature Tea assumed in the fit. The hori- 
zontal line represents the total ^-square value greater than 
the optimal value by unity, usually used as an error criterion. 
The corresponding plot of the total x-square values obtained 
in estimating Tsa is also presented in the inset. 
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FIG. 9. (Color online) The temperature and size depen- 
dence of the ratio of the CG and the SG correlation lengths 

tially an equal-weight sum of CG correlation functions, 
J cIxCk (x) . Since (x) in the large- a; region, which 
should govern the true CG correlation length, is much 
smaller in magnitude than that in the small- a; region, 
the latter contribution not playing an essential role in the 
CG correlation length might make a major contribution 
to ^iid mask or obscure the asymptotic be- 

havior of the CG correlation length. Thus, the standard 
OZ form may be inappropriate to treat the CG correla- 
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FIG. 10. (Color online) The spatial correlation functions 
of the SG and of the CG overlaps for several temperatures 
around Tea 0.308 and Tsa ^ 0.274. The upper curves 
are for the spin overlap and the lower ones are for the chiral 
one. The lattice is L = 16 under periodic BC. Total number 
of samples is Ns = 500. Error bars are omitted. Note the 
logarithmic scale of the ordinate. 

tion length, particularly when the system size is small. A 
similar observation was also made in the 4D Ising SG in 
a magnetic field^^. To overcome this difficulty, another 
dimensionless quantity was already proposed®*, but the 
examination of the quantity is beyond the scope of the 
present paper. 

Next, we turn to the quantities which probe the phase- 
space structure of the ordered state, including the overlap 
distribution and the non-self-averageness parameter. In 
Fig. [11] we show the distributions of the spin and the 
chiral overlaps at a temperature T = 0.2792, which lies 
below TcG and very close to (slightly above) Tsg- For 
L > 10, the chiral-overlap distribution shown in the 
right panel of Fig. [TT] exhibits a central peak in addition 
to the side peaks corresponding to the CG EA order pa- 
rameter ztq^Q. As the system size L increases, all the 
peaks grow in their height and become narrower in their 
width. This implies that all these peaks would remain 
in the thermodynamic limit. These features are noth- 
ing but the character of the IRSB, and are consistent 
with the occurrence of a negative dip in gcG- Similar 
behaviors were observed in several types of Heisenberg 
and XY SGs before^i^i^, but the side peaks observed 
in our present simulation for the SG seem sharper 
than those observed in the Heisenberg SG. 

By contrast, the spin-overlap distribution shown in the 
left panel of Fig.[TT]exhibits a shoulder-like structure only 
for small sizes, which tends to be suppressed as the sys- 
tem size increases. As the growing side peaks located at 
i'ZsG ^''^ expected in the SG ordered phase'*^, this obser- 
vation is consistent with the absence of the SG order at 
this temperature T — 0.2792, which is indeed compatible 
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FIG. 11. (Color online) The spin diagonal overlap distribution (left), and the chiral overlap distribution (right), at a temperature 
T = 0.2792 which lies in the CG ordered state, i.e. below Tea ^ 0.308 but slightly above Tsg ^ 0.274. A typical IRSB behavior 
is observed in the chiral overlap distribution. 



with our estimate above, Tsg = 0-274_q;q32. 

In Fig. [121 the SG and CG non-self-averagingness A 
parameters are plotted against the temperature for vari- 
ous system sizes. As can be seen from the figure, Aqg of 
different sizes intersect around T ~ Tea- A prominent 
peak is observed on the lower temperature side, which 
grows as the system size increases. This suggests that the 
self-averageness of the system is broken below Tcg and 
that the CG ordered phase is non-self-averaging. This is 
quite consistent with the IRSB nature of the CG phase, 
as already signaled by the central peak in and by the 
negative dip of gca- The parameter AgQ exhibits a be- 
havior similar to Aqg, accompanied with an intersection 
and a peak around T — Tea, even though the SG or- 
der is still absent at T = Tqc- As already noted in Sec. 
IIIH this is not surprising since a finite Asg just means 
the non-self-averageness of XSG which originates from the 
CG transition involving the phase-space narrowing asso- 
ciated with the IRSB. Hence, the intersection and the 
peak of Asg around Tcg is completely compatible with 
the spin-chirality decoupling. 

We also show the SG and CG G parameters in Fig.lT^ 
The observed behaviors are quite similar to the corre- 
sponding ones of the A parameters. Note that, although 
the intersections of the A and G parameters can also be 
used in estimating the transition temperature in princi- 
ple, the data tend to be noisy and not suited to precisely 
locate Tcg- The same suggestion was made for the 3D 
Ising SG^'S^ and for the 3D Heisenberg SG§5. 

V. CRITICAL PROPERTIES 

In this section, we study the critical properties of the 
SG and the CG transitions based on the finite-size scaling 
analysis. To estimate the CG critical exponents, we use 
both the CG Binder parameter and the CG susceptibility. 
On the other hand, we use the SG susceptibility only to 



estimate the SG critical exponents, the reason of which 
will be explained below. 

Let us start from the CG criticality. The standard 
finite-size scaling forms of the CG Binder parameter gcc 
and of the CG susceptibility XCG are given by 

gcG^x({T-TcG)L^''"'''). (24) 
Xcg = L^-''''^y({T-Tcg)L^'''''''). (25) 

where X and Y are appropriate scaling functions. 

For the CG susceptibility XCG, a good data collapse 
can be obtained by the two-parameter fits based on Eq. 

with vcG = 1-3 and -qcG = 0.23. Note that 
these values are obtained via the Baysian scaling anal- 
ysis (BSA), which enables us to estimate the critical ex- 
ponents in an unbiased way^, after fixing the transition 
temperature to the value obtained in the previous sec- 
tion Tcg = 0.308. By changing the assumed value of 
Tcg in the range of the associated error bar, we obtain 
the error bars of vcc and rjcc as vcg ~ 1.3 ± 0.1 and 
r]cG = 0.23 ±0.22. 

The crossing temperature of the CG Binder parameter 
gcG still exhibits appreciable size dependence, indicating 
the necessity to invoke the correlation-to-scaling term in 
the finite-size scaling. With the correction term, the scal- 
ing form of gcG is modified as 

gcG = A ((T - Tcg) L^/^^^^) (l + aL-"^«) , (26) 

where a is a numerical constant. The BSA analysis based 
on this form yields i^cG = 1-4 ± 0.1 and ujcg — 0.46 ± 
0.12, which is consistent with i/cc — 1-3 ± 0.1 estimated 
above from the CG order parameter. These values lead 
to (1/z^cg) + ^cg ~ 1-2, which is slightly larger than, 
but not incompatible with 9cg = 0-88 ± 0.03 obtained 
in the previous section. The resultant finite-size scaling 
plot is given in the left panel of Fig. [TH 
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FIG. 12. (Color online) The temperature and size dependence of the non-self-averageness A parameters of the SG (left), and 
of the CG (right). 
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FIG. 13. (Color online) The temperature and size dependence of the G parameters of the SG (left), and of the CG (right). 
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FIG. 14. (Color online) Finite-size-scaling plots of the CG Binder parameter (left), and of the CG susceptibility (right), where 
the correction-to-scaling effect is taken into account. The CG transition temperature is fixed to Tea ~ 0.308. Best fit is 
obtained for i/cG = 1-4 and ljcg = 0.46 for gcG (left), and ucc ~ 1-3, rjca = 0.29 and cjcg ~ 0.46 for XCG (right). 
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We examine the scaling form with the correction term 
also for the CG susceptibility, 

XCG = L^-^cGf (^(T - Tea) L^'^^^) (l + aL~^^^) .(27) 

Based on this form, we get vca = 1.3 ± 0.1 and rjcG = 
0.29 ± 0.12. These values are consistent with the val- 
ues obtained above without invoking the correction term. 
The resultant finite-size scaling plot is given in the right 
panel of Fig. [TH 

Next, we move to the SG criticality. As in the CG 
case, the standard scaling form of the SG susceptibility 
is given by 

XSG = L'-^'^Y ((T - Tsg) L^''^'^) . (28) 

A good data collapse of the SG susceptibility is ob- 
tained based on this form, with the resultant exponents 
vsG = l-23to TO and -qsG = -0.42l°;^^. The associated 
scaling plot is given in Fig. [151 Although we can obtain 



The estimated CG critical exponents are compatible 
with the ones reported in earlier literature on the 3D XY 
Sg46^^^^ They are also quite close to the values of 
the 3D Heisenberg SG^-, whereas they are largely differ- 
ent from the values of the 3D Ising SG, v — 2.5 ^ 2.7 and 
rj — —0.38 ^ — Q.40i2iii in spite of a common Z2 symme- 
try between the Ising spin and the chirality of the present 
model. The types of the RSB in the Ising SG and that 
in the XY or the Heisenberg SG may explain this differ- 
ence, i.e., the full RSB in the former versus the IRSB in 
the latter. The phase-space structure of the XFand the 
Heisenberg SGs is essentially different from that of the 
Ising SG. Such speculation also leads to another ques- 
tion: what causes the difference in the RSB types among 
the Ising, the XY and the Heisenberg SGs ? A possible 
explanation might be that the chirality-chirality interac- 
tion has a long-range nature different from the Ising one. 
Further study is needed to clarify these points. 
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FIG. 15. (Color online) Finite-size-scaling plot of the SG 
susceptibility without the correction-to-scaling term. The SG 
transition temperature is fixed to Tea = 0.274. Best fit is 
obtained for i^sG = 1-23 and rjsG ~ —0.42. 

an estimate of wsg by combining the value of vsG and 
that of 9sG — 1 -10-^55 estimated in the previous section, 
the result is inaccurate and poor, unfortunately. Other 
quantities such as the SG correlation length or the SG 
Binder parameter do not resolve this problem. Indeed, 
we also have examined the finite-size scaling of £,sg/L 
with the correction term, but ended up with an unphys- 
ical result of negative cusG- This inadequacy may partly 
be due to the fact that the estimated Tsg is located out 
of the range of the simulated temperature range. 

Summarizing the above results, we finally quote as our 
best estimates of the CG exponents as 

i^CG = 1-4 ± 0.1, 77CG = 0.29 ± 0.12, (29) 

and the SG exponents as 

i^SG = 1.23±°:j; vsG = -0.42^°:i2. (30) 



VI. SUMMARY AND DISCUSSIONS 

In this paper, we studied equilibrium ordering prop- 
erties of the 3D isotropic XY SG by means of exten- 
sive MC simulations, up to the linear size L = 40. Ex- 
amining various physical quantities including the glass 
order parameter, the Binder parameter, the correlation- 
length ratio and the overlap distribution function, we suc- 
ceeded in giving reasonable numerical evidence that the 
SG and the CG transitions occur at two different temper- 
atures. The estimated SG and CG critical temperatures 
are Tsg = 0.274j:J] °^^ ^^^-^^ j,^^ ^ o.308 ± 0.005, respec- 
tively, and the difference in the two transition tempera- 
tures is about 10 percent, which is comparable with the 
difference observed in the 3D Heisenberg SG. Our con- 
clusion of the occurrence of the spin-chirality decoupling 
in the model is in apparent contrast to that of the recent 
simulation by Pixley and Young^-. The main cause of the 
difference is that their analysis was mainly based on the 
correlation-length ratios. We have also confirmed that, as 
observed by Pixley and Young, the crossing points of the 
CG correlation-length ratio behave in not much different 
manner from the SG ones, which are seemingly consis- 
tent with a simultaneous SG and CG transition. Our 
present quantitative analysis, however, revealed that the 
extrapolated Tcross{L) of the ^cg/L lead to an unphysi- 
cal estimate of Tcgi *-e-, a negative one, whereas that of 
Csg/L lead to a reasonable estimate of Tsg- It turned 
out that such a behavior of the CG correlation-length ra- 
tio is inconsistent with that of other quantities such as 
the CG order parameter and the CG Binder parameter. 
In particular, the size dependence of the glass order pa- 
rameter speaks for the successive CG and SG transitions 
occurring at two different temperatures, as seen in Fig. 
[3l Our attitude is that the order parameter among var- 
ious quantities is expected to give a stable and reliable 
result, particularly when it is corroborated by the Binder 
parameter as we observed here, since it is the most fun- 
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damental quantity in describing phase transition. We 
then argued that the observed ih behavior of the finite- 
size CG correlation length, which is defined based on the 
OZ form, may be due to the special character of the CG 
spatial correlations deviating significantly from the stan- 
dard OZ form at short length and were susceptible to the 
pronounced finite-size effect. 

The critical properties of the SG and the CG orderings 
were also examined by means of the finite-size-scaling 
analysis. By controlling the correction-to-scaling effect, 
we obtained the CG critical exponents as vcg = 1-4 ± 
0.1 and rjcG — 0.29 ± 0.12. These values are close to 
the corresponding Heisenberg SG values, but are quite 
different from the Ising SG values in spite of a common Z2 
symmetry. The SG critical exponents were also estimated 
to be vsG = l-23lo;j!)g and rjsG = — 0.42to;27, which were 
consistent with the earlier estimates. 

The RSB nature of the ordered state was probed via 
the Binder parameter, the overlap distribution and the 
non-self-averageness parameter. All the quantities con- 
sistently point to the IRSB in the model. Physical signif- 
icance of the IRSB feature was already discussed in the 
Heisenberg case from several perspectives^. It would 
also be interesting to examine the possible IRSB proper- 
ties in real materials related to the XFSG such as granu- 
lar cuprate superconductors^^"— i^i^, together with the 
successive transitions and the associated critical proper- 
ties. 

As mentioned, the CG transition of the XY SG model 
and the Ising SG transition belong to different universal- 
ity classes in spite of a common Z2 symmetry between 
the chirality and the Ising spin variable. We speculate 
that this might originate from the difference in the type 
of the RSB in the two systems. The 3D Ising SG is be- 
lieved to exhibit a full RSB^i^>^, though some counter 
opinions also exist^"— . By contrast, the 3D A"y SG ex- 
hibits the IRSB. These different types of RSB may be 
related to the difference in the critical properties of the 
Ising SG transitions and the CG transition of the XY 
SG. Furthermore, since such a IRSB-like feature is also 
observed in the 3D Heisenberg SG^^, the above consid- 
eration suggests a possibility that the CG transitions of 
the XY and the Heisenberg SGs are actually the same, 
which are indeed consistent with our present numerical 
results. 
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Appendix A: The behavior of the spin-glass Binder 
parameter in the thermodynamic limit 

In this appendix, we examine how the SG Binder pa- 
rameter gsG in the thermodynamic limit behaves across 
TcG smd TsG in the occurrence of the spin-chirality de- 
coupling. 

In the CG phase realized at Tsg < T < Tcg, 
the average spin overlap (qa/s) vanishes as in the high- 
temperature paramagnetic phase. Meanwhile, the Z2 
spin-reflection symmetry is spontaneously broken in the 
CG phase, which implies that the determinant of the 
spin-overlap tensor, detqap, takes a nonzero value. De- 
noting this symmetry-breaking bias coming from the CG 
order by h{T), we may write the distribution of the spin- 
overlap tensor gap as 

P({g„4) oc e- ^ '^^) -'^m dot ^^^^ 

Since the standard SG order is absent in the CG phase, 
the average of any simple moment such as (q^p^ vanishes 
in the thermodynamic limit, whereas that of the deter- 
minant remains to be nonzero. Using this distribution, 
the SG Binder parameter is calculated as 

gsG = -Ih^iT). (A2) 

Hence, we find that the SG Binder parameter takes a 
negative value in the region Tsg < T < Tgg in spite 
of the absence of the long-range SG order. Although 
the detailed form of h{T) is unclear, h{T)^ would grow 
continuously and monotonically from zero when the tem- 
perature T is decreased across T — Tcg- It should be 
noted that, in the CG state of the 3D Heisenberg SG, 
gsG still remains to be zero in sharp contrast to the XY 
SG case95. 

In the SG ordered state reahzed at T < TsG: 9sg 
would take a different form. If there would be no RSB 
in the SG ordered state, gsG in the thermodynamic limit 
would jump to unity below Tsg- If there occurs the IRSB 
as in the present model, gsG would take a nontrivial value 
not equal to unity below Tsg, eventually tending to unity 
in the T — >■ limit. We show in Fig. [16] a schematic 
shape of gsG in the thermodynamic limit expected in 
the present model. 

Fig. [in] may enable us to extract further information 
about the transition temperatures from the data of the 
SG Binder parameter gsG for larger sizes. We see sev- 
eral characteristic temperatures in gsG for larger sizes, 
displayed in the left panel of Fig. U] For example, we see 
in the data of L = 40 two extrema and an inflection point 
in between. In the thermodynamic limit, the extremum 
at a higher temperature would converge to Tcg, while 
the one at a lower temperature to Tsc- Unfortunately, 
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the extrema are still faint and visible only for larger lat- 
tices, making a reliable estimate of Tqg and Tsg diffi- 
cult. The inflection point would be located somewhere 
between Tsg and Tcg in the thermodynamic limit, and 
thus, can be utilized to give a lower bound of Tcg and 
an upper bound of Tsg- Indeed, the inflection point of 
the L > 20 data lie around T ^ 0.31, which is compati- 
ble with our present estimates of Tsg = 0-2741q q32 and 
Tcg = 0.308 ±0.005. 



FIG. 16. A schematic form of the SG Binder parameter gsa 
in the thermodynamic limit, as expected in the present model. 
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